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Resonant tunneling in an open mesoscopic quantum dot is proposed as a vehicle to realize a 
tunable Fermi-edge resonance with variable coupling strength. We solve the x-ray edge problem 
for a generic nonseparable scatterer and apply it to describe tunneling in a quantum dot. The 
tunneling current power law exponent is linked to the S-matrix of the dot. The control of scattering 
by varying the dot shape and coupling to the leads allows to explore a wide range of exponents. 
Transport properties, such as weak localization, mesoscopic conductance fluctuations, and sensitivity 
to Wigner-Dyson ensemble type, have their replicas in the Fermi-edge singularity. 



Quantum dots host a number of interesting quantum 
transport phenomena, such as Coulomb blockade [1,2], 
Kondo effect [3] weak localization and universal conduc- 
tance fluctuations [4]. Electron scattering inside the dot 
as well as the dot-lead coupling can be controlled exter- 
nally by gates, which makes it possible to reveal phe- 
nomena of interest by varying system parameters. Such 
tunability has been exploited to demonstrate [5-8] new 
exotic varieties of Kondo effect. In this article we pro- 
pose to employ mesoscopic dots, in a similar controllable 
fashion, to study the Fermi-edge singularity (FES). 

FES is a fundamental manifestation of many-body 
physics taking place when an electron with energy just 
above the Fermi level tunnels into a metal, while leav- 
ing a localized hole behind. After tunneling, the elec- 
tron forms a quasiresonance due to interaction with the 
hole. This strongly affects the transition rate which is 
typically found to be a power law function of Mahan- 
Nozieres-deDominicis form A(e) oc (e — ep)~ a . Similar 
to the Kondo problem, the FES problem [9,10] is one of 
few many-body problems exactly solvable in the nonper- 
turbative regime of strong interaction. 

First discovered in the 60's in the context of x-ray ab- 
sorption in metals [9,10], the FES physics has found many 
other applications. In 1992, Matveev and Larkin [11] 
considered resonant tunneling and predicted a power law 
singularity, identical to FES, as a function of the reso- 
nance position relative to the Fermi level. In this case, 
the exponent a in the tunneling I — V characteristic is 
controlled by interaction of the tunneling electron and 
localized hole. The latter is system-specific, and depends 
on scattering phases and screening via Friedel sum rule. 

Below we generalize the theory [11] to describe reso- 
nant tunneling into an open quantum dot. Chaotic scat- 
tering in the dot, returns the tunneling electron many 
times to the hole, which enhances the FES singularity 
and makes it 'tunable', i.e. scattering-dependent. (In a 
noninteracting mesoscopic system [12], multiple returns 
to a resonant level are known to produce weak local- 
ization and conductance fluctuations.) While charging 
effects may interfere with resonant tunneling [13,14], in 
open dots one can ignore charge fluctuations and focus on 
the interplay of scattering and interaction with localized 



hole which forms FES. 

Manifestations of FES have been observed in tun- 
nel junctions [15] and in low temperature telegraph 
noise [16]. The role of scattering in quantum dots can 
be studied by resonant tunneling spectroscopy [17]. 

The canonical theory [10], based on separable scatterer 
model, is difficult to adapt to mesoscopic scattering. The 
crucial problem arises from noncommutativity of the S- 
matrices before and after electron release in the dot, ren- 
dering the separable model, along with the bosonization 
approach [18] used to handle it, irrelevant. Our approach 
builds on the Yamada and Yosida theory [19] of An- 
derson orthogonality for multichannel nonseparable scat- 
terer, recently advanced by Muzykantskii et al. [20,21], 
as well as on Matveev phase shift approach [22] to charge 
fluctuations in quantum dots. 

The theory presented below yields an exact relation of 
the FES exponent with the quantum dot S-matrix and 
reveals that the exponent structure is similar to that in 
the separable scatterer case. The orthogonality catas- 
trophe due to Fermi sea shakeup by switching of charge 
state at tunneling accounts only for one, negative part of 
the FES exponent, while the leading, positive part arises 
from interaction in the final state. Applying the result 
to the quantum dot problem, we find that by varying the 
dot scattering parameters the exponent a can be tuned 
to any value in the weak or strong coupling regime. 

Our results for open dots complement the work on the 
orthogonality catastrophe [23-25] and FES [26] in closed 
quantum dots which use the exact one particle states and 
energies to express the many-body overlap and transi- 
tion rate. The enhancement of orthogonality by disorder, 
discussed by Gefen et al. [25], has the same underlying 
physics as our backscattering-enhanced FES. 

The geometry of interest is pictured in Fig. 1 (a). We 
consider tunneling from a small dot which holds few elec- 
trons and has a large charging energy, into an open meso- 
scopic dot. The latter is characterized by a N x N S- 
matrix, where N is the number of channels connecting 
the dot and the leads. The interaction of electrons in the 
open dot with the localized hole in the small dot is de- 
scribed [22] by the backscattering phase 6 in the channel 
connecting the two dots (Fig. 1 (b)). 
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FIG. 1. a) An open quantum dot weakly coupled to a small 
closed dot which holds a localized electron that can tunnel 
into the open dot. b) Relation of the open dot S-matrix 
S and an auxiliary matrix S is illustrated. The latter de- 
scribes the dot with an extra open channel added to incorpo- 
rate backscattering on the small dot charge state. 

How does the dot S-matrix depend on backscattering? 
The answer is found most easily by considering an auxil- 
iary scattering problem with an additional channel which 
describes the point contact between the dots. This de- 
fines an extended S-matrix S of size (N + 1) x (N + 1). 
The physical matrix S can be linked with the auxil- 
iary matrix S by imposing the quasiperidicity relation 

Mtv+i = e ~ 2t5u N ,l +i on tnc m an d ou t components of the 
added channel, and eliminating these components from 

the scattering relation uf ut ^ = SijU^ n \ We obtain 

Sij = S ij + § ^g^i, i,j = l...N (1) 

with r = 5(jv+i)(at+i) the backscattering amplitude in 
the extended picture. One can verify that S*, defined by 
(1), is unitary provided that S is unitary. The relation 
between S and S is illustrated graphically in Fig. 1 (b). 

We emphasize that the parameters r and 8 which 
appear together in Eq.(l) and below describe different 
physics and, in particular, arise on different length scales. 
The phase shift S is a constant determined by the effects 
within a screening length from localized hole. In contrast, 
the quantity r, describing transport in the interior of the 
dot, is sensitive to the dot shape, and thus is tunable. 

The utility of the 5-dependent S-matrix (1) can be as- 
sessed by using the result for the orthogonality catastro- 
phe with nonseparable scatterer. In the latter problem, 
one is interested in the overlap of the many-body ground 
states with different S — So,i- Yamada and Yosida [19] 
derived a remarkably simple formula for this overlap, 

<l|0)=cxp(-^lniv), f = -L\tr In 2 S.S, 1 ] ( 2 ) 

with N the number of particles per scattering channel. 

For two different phases S, 5' in (1), the compound 
matrix R = S(S')S^ 1 (8) that appears in (2) is equal to 

= <% + (^- 1) m = e !^r r 1 (3) 

with Vi — — M 2 ) -1 / 2 , * = 1— W) the normalized 

column of S. Remarkably, R differs from a unit matrix 



only by a matrix of rank one. It means that, despite the 
complex dependence (1) on the scattering phase S that 
appears to affect the entire N x N matrix S, the orthog- 
onality problem is effectively a single channel-like. The 
overlap (0|1) is described by Eq.(2) with 

/3=-Imln(^) (4) 
A 2tt \U(5')J 1 ; 

Thus /? depends on transport in the dot solely via the 
backscattering amplitude r. Both the modulus of r and 
its phase, being functions of the dot shape and dot-lead 
transmissions, are under experimental control, and thus 
P can be tuned to any value in the interval < (5 < 1. 

We analyze the FES problem below for the scatterers 
(1), and find a relation between the FES and the orthog- 
onality exponents identical to the single channel case, 

a = 2(3-f3 2 (5) 

This is not entirely unexpected, given that the above 
analysis reveals hidden single-channel character of the or- 
thogonality problem. However, since the canonical FES 
theory [9,10] is limited to the separable scatterer situa- 
tion, the relation (5) cannot be deduced directly. Instead, 
we shall develop an approach for a generic S-matrix, and 
then specialize to the quantum dot case (1). 

Turning to the analysis of the FES problem, we con- 
sider the tunneling electron Green's function [9,10] 

G(r) = tr ($(0)e- iWlT $ + (0)e iWoT p e ) (6) 

with interaction included in the Hamiltonians Hi t 2 which 
describe electron scattering by the charged/uncharged lo- 
calized state. Here V>(t) = J2 a Ua ® a ( T ) * s the operator 
of a tunneling electron with a a labelled by energy and 
scattering channel, while p is electron density matrix 

Pe = 7F CX P (^-P^taa+a^j , ft- 1 = k B T (7) 

with the normalization factor Z = Y\ a (l + e - ' 56 "). The 
original approach [10] employs a diagrammatic expan- 
sion of (6) and, using the closed loop calculus, expresses 
it through electron Green's function describing time- 
dependent scattering at < t < t. Then the series for 
the Green's function are resummed in order to replace 
the scattering potential by the S-matrix. The resummed 
series, treated using Dyson equation in the time domain, 
lead to a singular integral equation that can be solved 
using a particular variety of the Wiener-Hopf method. 

Here we proceed differently, trying to avoid the dia- 
grammatic expansion altogether. This has a two-fold ad- 
vantage. Firstly, we shall be able to introduce the single 
particle S-matrices at an early stage of the calculation, 
thereby bypassing the rcsummation problem. Secondly, 
our approach will apply to noncommuting S-matrices. 

As a first step, we use the commutation relations 
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-ic a T iH T ~ + 



(8) 



a a Pe = e Pea a , K e 
to rewrite Eq. (6) as 

G(T) = 5>> a e^'e-^' r tr (e^ 1 T e m ° T Pe a+ a a ) (9) 



where a, a' label single particle band states. 



Next, we note that the quantities e 



-iHiT „iH T 



Pe 



are exponentials of operators quadratic in a a , dj, which 
allows to write their product as 



e- iHlT e m ° T p e = Z^exp 



where the operator w), defined by Eq. (10) and to be 
found in an explicit form below, acts in the single electron 
Hilbert space. With the help of the definition (10) the 
trace in Eq. (9) can be expressed through w as follows: 



tr (e- m ^e m ° T p e a+ a a ) = 



dct ( 1 



Z 



1 



r 1 , 

/ aa' 



which reduces the FES problem to analyzing the oper- 
ators 1 + e ±w . As we find shortly, the latter are re- 
lated to the single-particle S-matrix and energy distri- 
bution. The electron statistics is thus fully accounted for 
by the algebra involved in the construction of the opera- 
tor 1 + e w and its determinant, while the solution of the 
time-dependent scattering problem amounts to comput- 
ing the inverse (l + e~ u ') . The explicit separation of 
the many-body and the single-particle effects provides an 
efficient treatment of the FES problem. 

To make progress, we use the idea of Ref. [27] to link 
e w with single-particle quantities. From Baker-Hausdorff 
series for ln(e A e B ) in terms of multiple commutators of 
A and B, noting the correspondence between the com- 
mutator algebra of the many-body operators quadratic 
in a a , at, and the single-particle operators, we find 



e"' = e~ ihlT e ih ° T e~@* 



(11) 



Here the operators /io,i and e ^ e are related to the single 
particle Hamiltonian and density matrix (7) as 

H ,i = ^2(ho,i)aa>a+a a < , e _/3e = e~ Pec "5 aa > (12) 

aa' 

With the help of the result (11), defining n = (1 +e /3i )~ 1 , 
the determinant det (l + e 1 ") can be brought to the form 

det (1 + e*) - Zdet (l - n(e) + e - fhlT / h ° T n(e)) (13) 

The operator e~ thlT e lh ° T is represented most naturally in 
the basis of time-dependent scattering states constructed 
as wavepackets labeled by the time of arrival at the scat- 
terer. As Fig. 2 illustrates, the result of backward-and- 
forward time evolution is 



S = e ^ h ^e lhoT = S t t> 



R, < t < t 
1, else 



(14) 



with R = SiS Q 1 a compound S-matrix, and So,i the S- 
matrices for the charged/uncharged localized state. (In 
the single-channel case, R = e 2i ( 5 - 5 '\) Thus we rewrite 
Eq.(13) as det (l + e™) = Zdet (l + (S- l)nj . Simi- 

larly, the operator (1 + e 



'"' 1 1 becomes 



(l + e- 8 ) 1 = (n(e) + (l-n(e))S'- 1 ) ' n(e) (15) 



with n, S being operators in the Hilbert space of func- 
tions of time. 



SiSp \\fy^_: 



-ih,x 



S'o\ 



ih n x 



FIG. 2. Schematic forward and backward scattering time 
evolution e ~' h ^ T e th o r ; with ballistic transport outside the 
scattering region. The compound S-matrix R — SiSq 1 ac- 
counts for sequential scattering described by the S-matrices 
So,i corresponding to the hamiltonians ho t i. 

Thus the Green's function (9) is brought to the form 



G(t) = Le 



c <c 



det (1 + (S - l)n 



(16) 



L = J2e- ie ' T (Ue\(l-n)(n + S- 1 (l-n)) ' \u* e ,) (17) 



where u e = J2 a u aS(e — e a ) is a vector in channel space, 
and (...) includes summation over scattering channels. 
(A related determinant identity for e c has been known 
in the theory of counting statistics [28,20].) 

The factorization G(r) = Le c demonstrated for a gen- 
eral scattering problem with noncommuting So,i, pro- 
vides connection with Nozieres-deDominicis theory and 
generalizes it to nonseparable scatterers. The two factors 
in G(t), expressed through single-particle quantities, in 
the language of Ref. [10] correspond to the contributions 
of the open line and closed loop diagrams, respectively. 
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Backscattering phase, arg( r)/2n 

FIG. 3. The dependence (5), (4) of the tunneling exponent 
a on arg (r) for \r\ = 0, 0.1, ...0.9; Inset: The exponent a vs. 
\r\ for several values arg (r) marked by arrows (a — e). 



An explicit result for G(r) now follows by noting that, 
with respect to channel indices, the operators in (16), (17) 
are diagonal in the eigenbasis of R = SiSq 1 , where L is 
additive, while e c is multiplicative. Using the standard 
singular integral equation solution [10,19], we obtain 

3 

at t < h/T, with e 2m fy the eigenvalues of R, and £ ~ Ep 
the ultraviolet cutoff. The prefactor e _I<5M * describing the 
localized state energy offset can be discarded. 

In the case of our primary interest (3), R has only 
one nontrivial eigenvalue e 2m @ , given by (4), and Ui is 
an eigenvector Vi. This gives G(t) oc t'^~^\ leading 
directly to the result (5). The dependence of the FES 
exponent a on \r\ and argr is displayed in Fig. 3. 

The effect of mesoscopic fluctuations on FES can be 
described by drawing S from a Wigner-Dyson ensem- 
ble of matrices of size (N + 1) x (JV + 1), orthogonal, 
unitary or simplectic, depending on the symmetry. The 
backscattering amplitude r, being a diagonal matrix ele- 
ment of S, has a distribution [29] P(r) oc (1 — |r| 2 ) 7 with 
7 = N + 1, (N + 2)/2, 2N + 2 for the three ensembles. 
This generates an FES exponent distribution of width 
~ which is small at large N. For fixed modulus 

|r|, the change of the FES exponent can be of either sign 
depending on the phase 9 = argr (Fig. 3). The effect 
of scattering is particularly prominent at r approaching 
1, where the FES is strongly enhanced for the phase val- 
ues between 5 and 6', and suppressed otherwise, which 
corresponds to resonance formation in the dot. 

In summary, this work presents an exact solution of the 
Fermi-edge resonance problem for noncommuting scat- 
terers, relevant for tunneling in mesoscopic systems. We 
consider an application to resonant tunneling in open 
quantum dots and show that a resonance with tunable 



interaction strength, and thus with a variable power law 
exponent, can be realized. The resonance is strongly en- 
hanced by backscattering in a phase-sensitive fashion. 
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